Identification of dynamical Lie algebras for 
finite-level quantum control systems 



S. G. Schirmer, I. C. H. Pullen and A. I. Solomon 

Quantum Processes Group and Department of Applied Maths, The Open University, 
Walton Hall, Milton Keynes, MK7 6AA, United Kingdom 

E-mail: S . G . SchirmerOopen . ac . uk, I . C . H . PullenQopen . ac . uk, 
A.I. SolomonOopen .ac.uk 

Abstract. The problem of identifying the dynamical Lie algebras of finite-level 
quantum systems subject to external control is considered, with special emphasis 
on systems that are not completely controllable. In particular, it is shown that the 
dynamical Lie algebra for an N- level system with symmetrically coupled transitions, 
such as a system with equally spaced energy levels and uniform transition dipole 
moments, is a subalgebra of so(N) if N = 11 + 1, and a subalgebra of sp(£) if N = 2£. 
General criteria for obtaining either so(2£ + 1) or sp(£) are established. 
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1. Introduction 

In |]J we studied the problem of complete controllability of finite-level quantum systems 
with nearest-neighbour interactions. We showed that many quantum systems of physical 
interest are indeed completely controllable but that there are nevertheless systems with 
certain symmetries that are not completely controllable. This paper is devoted to 
identifying the dynamical Lie algebras for the latter systems. 

As in the previous paper, we consider the case of a driven quantum system, for 
which the interaction with the control field is linear, i.e., we assume that the total 
Hamiltonian of the system is 

H = H + f(t)H 1 , (1) 

where H is the internal system Hamiltonian and Hi represents the interaction of the 
system with the real control field /. We assume that Hq and Hi are Hermitian. For 
a finite-level quantum system there always exists a complete orthonormal set of energy 
eigenstates |n) such that H \n) = E n \n) and thus the internal Hamiltonian can be 
expanded in terms of the energy eigenfunctions |n), 

N N 

H = Y1 E n \n)(n\ = ^ E n e nn , (2) 

n=l n=l 

where e mn = \m)(n\ is an N x N matrix with elements (e mn )ki = o~ m kO~ni and E n are the 
energy levels of the system. The E n are real since H is Hermitian. We shall assume that 
the energy levels are ordered in a non- decreasing sequence, i.e., E\ E2 ^ - - - ^ En- 
Hence, the frequencies for transitions \n) —>■ \n + 1) are non-negative 

fi n = E n+1 - E n > 0, 1 < n < N - 1. (3) 

In the following it will be convenient to deal with trace-zero operators. Thus, if Hq has 
non- zero trace then we define the trace-zero operator 

H Q = H - [n^Tt(Ho)} I n , (4) 

which is equivalent to Hq up to addition of a constant multiple of the identity matrix 
In- Expanding the interaction Hamiltonian Hi with respect to the complete set of 
orthonormal energy eigenstates \n) leads to 

N 

H l = d m:n \m)(n\, 

m,n=l 

where the transition dipole moments rf m) „, which we assume real, satisfy d m ^ n = d n ^ m . In 
this paper we shall only be concerned with quantum systems for which the interaction 
with the control field is determined by transitions between adjacent energy levels, as 
is typical in the dipole approximation. It will also be assumed that there are no 'self- 
interactions', i.e., that the diagonal elements d n) n are zero for all n. Thus, letting 
d n = d n)n+ i for 1 < n < N — 1 we have 

JV-l N-l 

H i = d n(\ n )( n + l\ + \n+ l)(n|) = d n{e n ,n+i + e n +i, n )- (5) 

n=l n=l 
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2. Dynamical Lie algebras 

The operators iH^ and iH\ generate a Lie algebra C called the dynamical Lie algebra of 
the control system. This Lie algebra is important since it determines the Lie group S on 
which the control system evolves |2[]. Precisely speaking, the trajectories of the system 
subject to any control field are confined to the exponential image of the Lie algebra 
C. Knowledge of the dynamical Lie algebra thus enables us to determine the degree 
of controllability of a quantum system ||, f|, to identify reachable or non- reachable 
target states || |6|, and to determine whether a kinematical bound for an observable is 
dynamically accessible [[7], |f. 

The dynamical Lie algebra C generated by the operators iH and iifi defined in 
(0) and (H) is a real Lie algebra of N x N skew-Hermitian matrices, and the related Lie 
algebra £' generated by \H' Q and iH\ is a real Lie algebra of traceless, skew-Hermitian 
matrices. Thus, C is always a subalgebra of su(N). Since C is isomorphic to C @u{l) if 
Ti(H ) 7^ and C = CJ if Tr(fZo) = 0, it suffices to determine £'. It follows from classical 
results that a pair of skew-Hermitian matrices in su(N) almost always generates the full 
Lie algebra su(N) (see Lemma 4 in 0, for example). For the type of quantum systems 
considered in this paper, explicit criteria ensuring £' = su(N) have been established 

Theorem 1 Let do = d^ = and v m = 2d 2 n — d 2 mJrX — df n _ 1 for 1 < m < N — 1. The 
dynamical Lie algebra £' generated by iH' Q and iHi defined in (W and (j^) is su(N) if 
d m ^ 0, E m 7^ for 1 < m < N — 1, and one of the following criteria applies: 

(i) there exists \i v ^ such that fi m ^ fi p for m ^ p, or 

(ii) fi m = fi for 1 < m < N — 1 but there exists v p ^ such that v m ^ v p for m ^ p. 
If V = \N then d p -k ^ ±d p+ k for some k > is required as well. 

As has been shown in U, many quantum systems of physical interest indeed 
satisfy these criteria. However, there are systems of physical interest that do not meet 
these criteria. For instance, if any of the dipole moments d n vanish then the system 
decomposes into independent subsystems and its dynamical Lie algebra £' is a sum of 
subalgebras of su(N) ||. But even if all the d n are non- zero, the dynamical Lie algebra 
of the system may be a proper subalgebra of su(N), for example, if the transition 
frequencies /i n and the transition dipole moments d n satisfy 

\i n = fi N -n, d n = d N - n , 1 < Tl < N — 1, (6) 

as is the case for a system with iV equally spaced energy levels and uniform dipole 
moments. In the following we show that the dynamical Lie algebra £' of such a system 
is a subalgebra of so(2£ +1) if N = 2£ + 1, and a subalgebra of sp(£) if N = 2£, and 
give criteria ensuring £' = so(2£ + 1) or CJ — sp(£), respectively. In |Appendix 1J| , we 
also briefly discuss why the Lie algebra so{2£) does not arise for the systems considered 
in this paper. 
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3. The case N — 2£ + 1: dynamical Lie algebra so{2£ + 1) 
Consider a system with Hamiltonian H = Hq + f(t)Hi, where 

2£+l 21 

M = Y Enie n ,n, Mi = Y d n i(e n ^ n+ i + e n+ i >n ), (7) 

n=l n=l 

Ei < E 2 < ...E N , Ei 7^ E N , d n 7^ for all n, and the transition frequencies 
fi n = E n+ i — E n and transition dipole moments d n satisfy the symmetry relation (g). 
We shall prove that the Lie algebra £' is a subalgebra of so(2£ + 1), which is in general 
isomorphic to so{2£ + 1). 



3.1. C C so{2£ + 1) 

We show first that £' C so(2£ + 1). Let y n>m = i(e H)m + e mi „). Using d n = c^+i-n and 
2/m,n = J/n,m we can simplify 

i 

Mi = Y de+i-n(ye+2 —n,£+l—n + ye+ n ,£+n+i)- 

n=l 

To compute Hq, we note that E n = E\ + Z)"=i /V Thus, using [L n = fj, 2 e+i-n leads to 
Tr(F ) = (2£ + l)Ei + (2£ + 1) ]T 

8=1 

Hence, the energy levels £^ of H' Q are = and 

1 e+n e 

E'e + i- n = — Y E' e+ i +n = Y = Y V s 

s=£+l—n s=e+l s=£+l—n 

for 1 < n < £. Consequently, we have 

1 ( 1 V 

M' = X/ I ~ A*s I K e £+l-n,£+l-n — Q+l+n,£+l+ra)- 

n=l \ 8=i+l— n / 

Let cr be an isomorphism of the Hilbert space of pure states defined by 

/ K + 2-n), l<n<£ + l 
ff(|n>) = \ £ + 2<n<2£ + l (8) 

and set |m) = c(|n)) as well as i£ m = — Yf s =i+i-m A*a an d = ^+i-m- Then the 
representations of iifp and iifi with respect to the new basis \m) are 



m=l \ s=£+l— m / m=l 



(9) 



m=2 m=l 

with /i m and ?/ m as defined in ( |B.1| ) and (|B.3| ), respectively. Hence, M Q and 1H1 are both 
in so(2£ + 1) and thus the Lie algebra CJ they generate must be contained in so(2£+ 1). 

Since iH' Q and iifi in (|Sp contain a complete set of generators h m and y m for so(2£+l) 
(see |Appendix B|) , it is natural to expect that they generate the full Lie algebra so(2£+l). 
We shall prove that this is usually, but not inevitably, true. 
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Example 1 Consider a system of type (ffi) for £ = 3. If E m = m, 1 < m < 7 and 
d\ = c?6 = v^3? c?2 = <^5 = v^5 and d% = d^ = v^6 i/ien the basis change (Qj /ea<is to 

iH'q = — hi — 2/i2 — 3/i3 

ii/i = v^yi + y/by 2 + V3y3 

with h m and y m as defined in ( \B. 1\ ) and ( \B. 3\ ), respectively. Therefore, the Lie algebra 
£' generated by iH' and 1H1 is a subalgebra of so{7). However, it is easy to verify that 
£' 9^ so(7). Indeed, in this particular case £' is a three-dimensional subalgebra of so(7) 
spanned by iH , 1H1 and [iH , iHi] = y/6x\ + \fh%2 + V^s- 

Thus, for certain choices of the parameters E m and d m , the Lie algebra £' is a proper 
subalgebra of so{2£ + 1). 

3.2. Criteria for £' = so{2£ + 1) 

To find criteria that ensure £' = so(2£ + 1), consider the generic system 

t i 
lH o = H e m /i m , 1H1 = SmUm, e m 7^ 0, 5 m ^ Vm (10) 

m=l m=l 

with h m and y m as defined in (|B.1|) and ( |B.3|) , respectively. As before, iH and iHi are 
in so{2£ + 1) and hence the Lie algebra £' they generate must be contained in so(2£ + l). 

Theorem 2 Let u m = e m+ i — e m for 1 < m < £ and ujq = ei- The dynamical Lie 
algebra £' generated by the system H = H' Q + f(t)Hi with iH'q and iHi as in fiTty) is 
so(2£ + 1) ifuj 2 m ^ul for\<m<£. 

PROOF: Using the properties of the generators h m and y m leads to: 

V® =[[iH' ,iH 1 )M 

1 

m=l 

V« =[[rHiVt%iH> Q ]-u>UV {G) 

e-i 

m=l 

e-2 



m=l 



i-\ 

m=l 

By hypothesis 7^ for m > and 5i 7^ 0, u; = e x 7^ 0. Hence, all the factors in 
the last expression above are non-zero, i.e., we have y\ G £' and thus £' = so(2£ + 1) 
by lemma [l] of [Appendix B| . □ 
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If u;^ = Uq for some m > then a slight modification of the proof above leads to a 
residual term 

e 

y(o) = ^ 5 m y m = bmy m (11) 

m£M m=l 

where M. = {m : 1 < m < £, ui^-i = ^o) with S m = 5 m for m G M. and 5 m = 
otherwise. If the energy levels are either positive and ordered in a non-decreasing 
sequence, i.e., < e m < e m+ %, or negative and ordered in a non- increasing sequence, 
i.e., > e m > e m+ i, then = implies £J m _i = u for all m G M.. We shall only 

consider this case in the following. 

Theorem 3 Let v m = 25^ — S^+i — for 1 < m < I, where 5 = 5i, 5e + i = 0. The 
dynamical Lie algebra C generated by the system H = H' Q + f(t)Hi with iH' and iHi 
as in fcu\ ) is so{2£ + 1) if ' 0J m -\ = but v m ^ V\ for all m G M. — {1}. 

PROOF: Since io m -\ = for all m G Ai, we have 

Z =2-^\Y^]=j:C5 2 m+1 -~5 2 Jh m . 

m=l 

Suppose M. — {1} has I' elements labeled mi, 777-2 up to my. If v m 7^ v\ for all 
m G M — {1} then 

= 5i(t>i - u TOi ,)yi - 5 mk (v mk - v me ,)y mk 

= - U m ^,)Xi - ^ ^m fe (fm fc - V me ,)x mk 

^E[Z,lW]-V-/ (1) 

€'-2 _ 

= 5i(7Ji - U m< ,)(Vi - V^^yi - Y bm k { V m k ~ ~ V rr H t ^Vm k 

xW = [YU,Z]-v mel XM 



Sl(Vl - V mil )(V\ - V mel l )Xi - $m k {Vm k ~ %,)(% ~ "m,,^)^ 

k=l 



k=l 



shows that y% G C and hence C = so{2£ + 1) by lemma [l] of |Appendix B| . □ 



A similar argument shows that if there exists k G M such that v m 7^ Vk for 
m G M. U {1} but m y£ k, then y fc G Using the fact that the generators y m of 
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so(2£ + 1) are not diagonal with respect to the Cartan elements (see [Appendix lj| ), it 



generally follows that C contains all the generators x m and y m and thus C = so(2£ + 1) 
as well. An important special case of this type is a system with N = 2£ + 1 equally 
spaced energy levels and uniform transition dipole moments: 

Theorem 4 The dynamical Lie algebra CJ generated by the system H = H' Q + f{t)H\ 
with N = 2£ + 1 equally spaced energy levels u m = e± and uniform dipole moments 
S m = 5 is so(2£+ 1). 

Proof: Let = S'H^ and X« = e^xH'^Y^). Then h t = -2- 1 [X«, Y"«], 
yt = [h £ ,X^\ and x t = [y e , h £ \. Thus, x t , ye G C . Next, set = F« - y e+1 _ k and 

X (k+i) = X {k) _ X£+i _ k f or i < A; < £, and note that h e _ k = 2- l [Y ( - k+1 \ X^ k+ % 
ye-k = [h>t-k,X ( k+1 >] and x^k = [ye-k, he-k]- This shows that CJ contains all the 
generators x m and y m of so{2£ + 1). Hence, CJ = so{2£ + 1). □ 

3. 3. Application of the criteria 

We shall now return to the original system (0). Since we have assumed that the energy 
levels are ordered in a non-decreasing sequence, we have \i m > for all m and hence 

e 

E m = - - °> E ™ - E ™+i 1 <m<£ 

s=e+l—m 

i.e., the energy levels E m are negative and form a decreasing (non-increasing) sequence. 
Noting that e m = E m , we thus have u = e% = E\ = —jie+i < and 

i i 

w m = e m+l — e m = E m+ i — E m = ^ fl s — ^ fl s = —fl£- m +l — 0. 

s=£— m+1 l—m 

If fii ^ and n m ^ fig for m < £ then £ = so(2£ + 1) according to theorem 
independent of the choice of the dipole moments, provided that they are non-zero. 

If fi m = fi£ for some m < £ then theorem |3| applies and the dipole moments 
determine whether we have CJ = so(2£ + 1) or a proper subalgebra. In particular, 
the Lie algebra is so{2£ + 1) if the energy levels are equally spaced, \i m = ^ for all m, 
and the dipole moments d m are such that v m ^ v± for all m > 1, where 

„ _ 2d 2 -d 2 -d 2 - 2d 2 -d 2 -d 2 

for 1 < m < £ and v\ = d\ — d\ = d\ — dj_ l (with do = 0). For instance, 
setting d m = y/m for 1 < m < £ works since it gives V\ = £—{£—!) = 1 but 
v m = 2{£ + 1 — m) — {£ + 2 — m) — {£ — m) = for m > 1. Setting d m = 1 also gives 
CJ = so(2£ + 1) due to theorem |. 

However, recall that there are systems of type (^) whose Lie algebra £ is a 
proper subalgebra of so{2£ +1), as example |I| above clearly shows. Note that in this 
example, the energy levels are equally spaced and the dipole moments d m are such 
that vt = 2d§ - d\ = 6 - 5 = 1, v 2 = 2d\ - d\ - d§ = 2 x 5 - 6 - 3 = 1 and 
t>3 = 2d\ — d 2 , = 2 x 3 — 5 = 1, i.e., all the v m are equal and none of the theorems above 
are applicable. 
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4. The case N = 2£: dynamical Lie algebra sp(£) 

Let N = 2£ and consider the system H = Hq + f(t)Hi with 

21 21-1 

iH = E n ie n>n , lH^ = ^ d n y n>n+1 , (12) 

n=l n=l 

where fi n = and d n = d<it- n for I < n < £ — I. Note the symmetry of the system. 

Every transition frequency except fie occurs in pairs, and although fi£ may be different 
from all the other fi n , theorem [l] does not apply since N = 2£ and there is no k such 
that de-k ^ ±G?£+fc. In fact, we shall prove that the Lie algebra generated by iH' Q and 
iH\ is a subalgebra of sp(£), which is in general isomorphic to sp(£). 



4.1. C! c sp(£) 

To show that £ C sp(£), we note that c^-n = d n implies 

l-\ 

iHi = diyi t i + i + 22 d n (y n , n +i + y2£-n,2e+i-n)- 

n=l 

Furthermore, E n — E\ + S"=i anc ^ ^u-n — A*n f° r 1 < w < ^ — 1 leads to 

€-1 

Tr(# ) = (2^)^i + (2£) X; ^ + ^ 



Hence, Tr(if ) / (2^) = -Ei + Z)«=i + |w an d the energy levels of H' are £7^ = — §/^, 
E' l+l = \l±i and 

£-1 ,. e+n ,, £-1 

^£-n = - 2^ - Y ' = 2^ ^ s + y = + y • 

s=£— n ~ s=l+l ~ s=l—n ~ 

for 1 < n < I — 1. Thus, we have 

i-^CI = XZ — ( / Xs ) K e n,n — £2l+l~n,2l+l-n) 

n=l V ^ / 

Let a be an isomorphism of the Hilbert space of pure states defined by 

|n) 1 < n < 

-l) n - £ ~ 1 \3£+ 1 -n) £+l<n<2£ 



< J (\ n )) = \ < ^„-/-l,o, , 1 „ \ , ,i^„^o, ( 13 ) 



and set |m) = cr(|n)) as well as _E m = — |//^ — YslZm^s for 1 < m < £ — 1, ^ = — |/^, 
^ = dg. Then iifg and iifi have the following representations with respect to the new 
basis \m) 

1 (ill l ~ x \ e 

iH' Q = X — I ~7T + X ^ s J K e m,m — €-m+e,m+e) = X E m h m 

m=l \ s=m ) m=l I'l/H 

m=l m=l 

where h m and y m are as defined in flU.lfl and ( |(J.3| ), respectively, and we note that 
y m ,m+i = y m +i,m, y m+ e,m+e+i = y m +i+i,m+e- Hence, the dynamical Lie algebra generated 
by iH' and iH 1 must be a subalgebra of sp(£). 
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Since iHq and iHi in (|14|) contain a complete set of generators h m and y m for sp(£), 
it is natural to expect that they generate the full Lie algebra sp(£). We shall prove that 
this is true in most cases. However, as in case of so(2£ + 1), a proper subalgebra may 
also be generated. 

Example 2 Consider a system of type fiLl\ ) for £ = 3. // E m = m, 1 < m < 6 and 
d\ = d$ = \JT>, d 2 = d± = 2\/2 and d 3 = 3 then the basis change ftTSj ) leads to 

iH' Q = - 2.5/ii - 1.5h 2 - 0.5h 3 

Lffi = Vhy x + 2y/2y 2 + 3y 3 

with h m and y m as defined in (\C. 1\ ) and (\C. respectively. Therefore, the Lie algebra 
£ generated by iH' Q and iHi is a subalgebra of sp{3). However, it is easy to verify that 
C ^ sp(3). Indeed, CJ is a three-dimensional subalgebra of sp(3) spanned by iHq, iH\ 
and [iH , iifx] = — (v^^i + 2\[2x 2 + 3x 3 ). 

Thus, for certain choices of the parameters E m and d m , the Lie algebra C is a proper 
subalgebra of sp(£). 



4-2. Criteria for £ = sp{£) 

To find conditions that ensure CJ = sp(£), we consider the generic system 

t i 
\H' = £ m h m , \Hi = Y Smy m , e m ^0, 5 m ^ Vm (15) 

m=l m=l 

with h m and y m as in ( |C.1| ) and ( |C.3| ). Clearly, iH' Q and iHi are in sp(£). Hence the Lie 
algebra CJ they generate must be contained in sp(£). 

Theorem 5 Let u m = e m+ \ — e m for 1 < m < £ and u>e = 2e^. The dynamical Lie 
algebra CJ generated by iH' and iH 1 as in ( TH| ) is sp(£) if ' ^ ujJ for m < £. 



PROOF: Using the properties of the generators h m and y m leads to: 

t 

= J2 *m(^ - ul)y m 

m=2 

=[[iKv(%iH' ]-u;lvU 

i 



m=3 



m=l 

By hypothesis u;^ ^ ujJ for m < £ and 5i ^ 0. Hence, all the factors in the last expression 
above are non-zero, i.e., we have ye £ CJ and thus CJ = sp(£) by lemma @ of [Appendix 
0. □ 
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If = ujf for some m < £ then a modification of the proof above leads to a residual 

term 

i 

Y^ ^ = ^ ] ^rnVm ^ ] 3m,ym 
m£M m=l 

where M. = {m : 1 < m < £, — ujf } and S m = 5 m for m G M. and 5 m = otherwise. 

If the energy levels e are negative and ordered in an increasing (non-decreasing) 
sequence then = uj implies u m = —u e = — 2q > for all m G M.. We shall only 
consider this case in the following. 

Theorem 6 Let v m = 1b 2 m — <5„ +1 — for 1 < m < £ and 6t+i = <^-i, 5o = 0. The 
dynamical Lie algebra C generated by the system H = H' Q + f(t)H\ with iH'q and iH\ 
as in ( |75| j is sp(£) if /i m = —fAi = —2e e but v m ^ v% for all m G M. — {£}. 

Proof: Let X<® = -fiJ^iHo, Y^] and 

m=l 

Suppose M. — {£} has £' elements labeled mi, rri2 up to me> and let m^+i = £. Then 
= [Z,X®] - v mi Y<® = £ ~5 mk (v mk - v mi )y mk 

k=2 

e'+i _ 

X (D = [y (o) ; z] _ V ^ X {0) = ~5 mk (v mk - v mi )x mk 

k=2 
I'+l _ 

F {2) = [Z,X {1) ] - V m2 Y W = S m k ( V m k ~ V)K - V m 2 )ym k 

k=3 

X^ = \Y^,Z] - v m2 X^ = - v mi )(v, 

k=3 



Y^ = f[6 t (v e -v mk )y t 
fc=i 



m k V rrl 2)X rnk 



shows that yi G CJ and hence C = sp(£) by lemma of |Appendix (J| . □ 



A similar argument shows that if there exists k G M such that v m ^ Vk for m,k G M. 
but m k, then y*. G Using the fact that the generators y m of sp(£) are not diagonal 
with respect to the Cartan elements (see [Appendix C| ), it generally follows that C 



contains all the generators x m and y m and thus C = sp{£) as well. An important 
special case of this type is a system with N = 2£ equally spaced energy levels and 
uniform transition dipole moments: 

Theorem 7 The dynamical Lie algebra C generated by a system H = H' + f(t)Hi with 
N = 2£ equally spaced energy levels and uniform dipole moments is sp{£). 
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PROOF: We have xt = uo y e = —2 1 [iH' ,x e \ and h £ = —2 l [x^yA. 

Next, set = Y^-^ - y e+1 _ k and Z&> = Z^-^ - h e+1 _ k for 1 < k < £, with 
y(°) = and Z(°3 = \H' Q , and note that x t - k = [Z^,Y^\, y^ k = -2~ X [Z^ , x e . k ] 

and hi_ k = —2~ x \xi^ k ,yi^ k \. This shows that £ contains all the generators x m , y m of 
sp(£). Hence, CJ = sp(£). □ 



4-3. Application of the criteria 



Let us now return to the original system (0). Since we have assumed that the energy 



levels E m of the system are ordered in a non-decreasing sequence, i.e., E m < E m+ i for 
all m, we have \i m > 0, and hence 



E m = ~ + J2 l^s) < 0, E m < E m+1 , 1 < m < 




i.e., the energy levels E m are negative and form an increasing sequence. Noting that 
E m = e m we thus have u) m = e m+1 — e m = \i m >0forl<m<£ — 1 and u>i = 2ei = 

Thus, if Hi ^ and /i m 7^ /Z£ for m < £ then £' = according to theorem |, 

independent of the choice of the dipole moments, provided that they are non-zero. If 
— He for some m < £ then theorem |6] applies and the dipole moments determine 
whether we have C = sp(£) or a proper subalgebra. 

In particular, the Lie algebra is sp(£) if the energy levels are equally spaced, fi m = fi£ 
for all m, and the dipole moments d m 7^ 0, 1 < m < £, are such that there exists k, 
1 < k < £, so that v m 7^ v k for all m 7^ k, where 

v _ 9J 2 — A 2 _ J 2 — 9r/ 2 — r/ 2 — r/ 2 1 < m < / 

(with <i = 0) and ve = 2d} — 2dj_ 1 . 

For instance, setting d m = i/m for 1 < m < £ works since it gives v m = 
2m — (m — 1) — (m + 1) = for 1 < m < £ but = 2. Similarly, setting d m — 1 
for 1 < m < £ also gives CJ = sp(£) due to theorem [7|. 

However, recall that there are systems of type ([T^) whose Lie algebra CJ is a 
proper subalgebra of sp(£), as example ^ above clearly shows. Note that in this 
example, the energy levels are equally spaced and the dipole moments d m are such 
that Vi = 2d\ - d\ = 2 x 5 - 8 = 2, v 2 = 2d\ - d\ - d\ = 2 x 8 - 5 - 9 = 2 and 
t>3 = 2d\ — 2d\ = 2(9 — 8) = 2, i.e., all the v m are equal and none of the theorems above 
are applicable. 



5. Conclusion 



Our analysis of finite-dimensional, non-decomposable driven quantum systems with 
nearest neighbour interactions and non-zero dipole moments shows that the dynamical 
Lie algebra CJ generated by the trace-zero part of the internal Hamiltonian and the 
interaction Hamiltonian of the system is either su(N), so(N), sp(^N), or a simple 
subalgebra of these. Although by far the most common case is su(N), which corresponds 
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to density matrix / observable controllability and usually complete controllability ||, 
certain symmetries of the controlled transitions can destroy or reduce the controllability 
of the system. 

Precisely, we showed that the dynamical Lie algebra C of a system with 
symmetrically coupled transitions is a subalgebra of so(2£ + 1) if the system has an 
odd number of energy levels (where degenerate levels are to be counted according to 
multiplicity) and a subalgebra of sp(£) if the system has an even number of energy levels. 
Moreover, we established criteria which guarantee in most cases that the dynamical Lie 
algebra is actually isomorphic to either so(2£ + 1) or sp(£). In particular, the dynamical 
Lie algebra of a system with equally spaced energy levels and uniform transition dipole 
moments is so(N) if N = 2£ + 1, and sp (|JV) UN— 2£. 

Despite the rather technical nature of the results presented in this paper, we would 
like to emphasize that the identification of the dynamical Lie algebra is a crucial first 
step towards identification of reachable and non-reachable target states for systems 
that are not essentially controllable. Furthermore, knowledge about the structure of the 
dynamical Lie algebra can be used to develop efficient control schemes for these systems. 



Appendix A. The Lie algebra su(N) 

A standard basis representation for the Lie algebra su(N) in terms of trace- zero, skew- 
Hermitian N x N matrices is (see, for example [[K]]) 

Vm,n = i(^m,n &n,rn)i (A • ^ ) 

hm = l\&m,m Sm4-l,m+l)) 

where 1 < m < N—l, m < n < N and i = \J — 1. There are £ = N — 1 generators h m and 
\£{£+ 1) generators of type x m>n and y m>n each. Hence, the total number of generators 
is N 2 — 1 and thus the dimension of the Lie algebra su(N) is iV 2 — 1. A nice discussion 
of controllability of iV-level quantum systems in terms of root space decompositions of 
su(N) can be found in |§. 



Appendix B. The Lie algebra so{2£ + 1) 

so(N) usually refers to the real Lie algebra of trace-zero, anti-symmetric matrices (lCj. 
However, since we are dealing with subalgebras of su(N) generated by iV x N skew- 
Hermitian matrices, we require a representation of so(N) in terms of trace-zero, skew- 
Hermitian matrices. For N — 2£ + 1, the standard representation of the complex Lie 



algebra Bi |T1J leads to the following skew- Hermit ian basis for the real Lie algebra 
so(2£ + 1): 
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h-m = i( e m+l,m+l — e m+£+l,m+£+l) 

X €m — Xl,m+1 ~ Xm+£+l,l 

Ue m = Vl,m+1 ~ Um+e+1,1 

X e m +t n = x m+£+l,n+l ~ x n+l+l,m+l (B-l) 

Ve m +en = y m +l+l,n+l ~ 2/n+£+l,m+l 

X e m -e n = x n+l,m+l ~ x m+l+l,n+£+l 

2/e m -e„ = Un+l,m+l — Vm+e+l,n+l+l 

where 1 < m < £ and m < n < £. Since there are £ elements h m , x tm and y em each, as 
well as - 1) elements x em+£n , y em +e n , x e m -e n and y 6m -e n each, the total number of 
basis elements is £{2£ + 1). Thus, the dimension of so{2£ + 1) is £{2£ + 1). Using the 
general commutation relations 



(B.2) 



(B.3) 



[ X e m l 2/e,J Htm-tri 2/e OT - < „ 

k, 2/eJ = -2/im, [Xe m ±e n ,ye m ±e„] = ~^{h m ± K), 

[h m , X tjn ± tn ] = —y ern ±e n , [h m , 2/e m ±e n ] = x e m ±e n 

for m n, shows that the elements x m and y m with 

x l X eii x m+l x e m — e m j r ii 1 — ^1 — £ 1) 

yi = , 2/m+l = ye m -e m+1 , 1 < m < £ - 1, 

are not diagonal with respect to the Cartan elements h m of the Lie algebra and generate 
the full Lie algebra so(2£ + 1). Furthermore, it generally suffices to prove that the 
Lie algebra CJ generated by iH' and iH\ as in ( |i0|) contains one of these elements to 
conclude that CJ = so(2£+ 1). We shall demonstrate this explicitly for the case y\ G L' . 

Lemma 1 Let CJ be the Lie algebra generated by iH' Q and iHi as defined in pTL 1 . If 
2/i G £ then x m , y m G CJ for 1 < m < £ and hence CJ = so{2£ + 1). 

PROOF: Using ( B.2|) shows that y 1 G C implies [iHo>2/i] = e i x i an d [^i, 2/i] = 2h x \ thus 
x\,h\ G CJ . Furthermore, we have 



-I in 



m=2 

I 

m=2 

I 

= -[iH , iili] + eiSiX! = XX e ™ _ e m -J)b~ m x m 

m=2 

I 

[Z {1 \ Y^] = -e 2 5 2 x 2 - J2(e m - e m _i)5 m 2 r . 



m=3 



which shows that + [Z {1 \Y^\ = -e x b 2 x 2 , i.e., x 2 G C, and [Z ( - 1 \x 2 ] = e 2 y 2 , 
[#2,2/2] = 2(h 2 — h\) implies y 2l h 2 G CJ . In general, defining recursively 

Z (k) = Z (fc-1) _ €khki yik) = y(k-l) _ X (k) = X (k-1) _ (efe _ €k l)SkXk 

shows that X (fe) + = -e k S k+1 x k +i, [Z (fc) , x k+1 ] = e k+1 y k+1 and [x k+ i,y k+1 ] = 

2(/ife+i — h k ). Thus, x k+ i, y k+ i and h k+ i are in CJ for k — 2, 3, ...,£ — 1. □ 
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Appendix C. The Lie algebra sp(£) 



A basis representation for the Lie algebra sp(£) for iV = 2£ in terms of trace-zero, 



skew-Hermitian N x N matrices can be derived from the standard basis for CV |TT 

hm i(^m,m &m+l;m+t) 

X2t m •^m+£,m 

y2t m = ym+e,m 

^e m +e n •^m+t,n •^n+£,m (0.1 J 

Ve m +e n — Vm+£,n + Vn-\-l,m 
•^e m —e n •En,m %m+l,n-\-l 

Hem— en = Vn,m ym+£,n+£, 

where 1 < m < I and m < n < I. Since there are I elements h m , x 2tm and y2e m each, 
as well as \£{£ - 1) elements x £m+e „, y £m + £n , x Em _ en and y em -e n each, the total number 
of basis elements is £(2£ + 1) and the dimension of sp(£) is thus £(2£ + 1). Using the 
general commutation relations 

[ X 2e n i %e m — £„] ^em+En) \^2e n i Hem — e„ ] Df-m+tni 

Ffm+enl^m-en] 2(x2e m X2t n ), [^tm+en! ?/em~£n] 2(?/2e m + Z/2e„)j 2^ 
[^2e m ,2/2e m ] = ~2/i m , [^6 m ±6 ra , 2/e m ±eJ = -2(/l m ± K), 

[h m , X em ± €n ] — —ye m ±e n -, [hm,y €m ± €n \ ~ Xe m ±e n 

for m ^ n, shows that the elements x m and y m with 

x m = x em - tm+1 , l<m<£-l, x £ = x 2ev 
y m = Ve m -e m+1 , l<m<£-l, y e = y 2et . 

are not diagonal with respect to the Cartan elements h m of the Lie algebra and generate 
the full Lie algebra sp(£). Again, it therefore generally suffices to prove that the Lie 
algebra £ generated by iH'q and iHi as in (|1|) contains one of these elements to conclude 
that £ = sp(£). We demonstrate this explicitly for the case ye G V . 

Lemma 2 Let CJ be the Lie algebra generated by iH' Q and iHi defined in (|73[). Ifyt G CJ 
then x m , y m G CJ for 1 < m < £, hence £ = sp(£). 

PROOF: Using (|C.2| ) shows that yi G CJ implies [iH^y^] = 2eiXi and \xi,yj\ = 2hf, thus 
xg, he G CJ . Furthermore, we have 

l-\ 

' z m IL m 



m=l 

£-1 

= iifi - 6ty t = 8 mVm 

m=l 
£-1 

X (1) = -[iH' Q , iHt] + 2e £ 8 £ x £ = ^ (Wi ~ ^m)S m x m 

m=l 

£-2 



m=l 
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which shows that + [,Z*W,YW] = e^_ia^_i, [^W,^] = — ee-iye-i, and 

[a;^_i,2/£_i] = 2(hi — thus a^_i,?/£_i, G In general, defining recursively 

X (*+D = X (k) _ (Q _ fc+i _ ei _ k ) St _ k x e _ k 

shows that x( fc+1 ) + [Z( fc+1 ),y( fe+1 )] = e t - k 5t-. k - 1 x i -. k - 1 , [Z^^x,^] = -e^m^ 
and [x^_ fc _i,^_ fe _i] = 2(h £ _ k -h i _ k _ 1 ). Thus, we have indeed a^_jt_i, and /^_ fc _i 

in C for Jfe = 1,2,... ,£-2. □ 



Appendix D. The Lie algebra so(2£) 



Using the standard representation for the complex Lie algebra D £ p| , we can derive 
the following skew-Hermitian basis for so(2£): 

Vem+en = Urn+e,n ~ Vn+l,m (D-l) 
Ue m —e n 2/n,m ym+t,n+£ 

where 1 < m < £ and m < n < £. There are £ elements h m , as well as \£{£ — 1) 
elements x em+en , y em+en , x tm - en and y tm -e n each, i.e., the total number of basis elements 
is £{2£ - 1). Thus, the dimension of so{2£) is £{2£ - 1). 

To see why there is no (2£) -level system with H = Hq + f(t)H\, where iHo and iH\ 
are as defined in (Q) and (|5]), respectively, such that £J = so(2£), note that 

x m = x e m -e m+ ii 1 < m < £ — 1, X£ = X ei _ 1+ei , ^ 

y m = y em - em+1 , l<m<£-l, y e = y £e _ 1+£e . 

forms a minimal, complete set of generators for so(2£) if £ > 2. Each of the £ generators 
y m has four distinct, non-zero entries, which corresponds to a total of 4£ non-zero entries. 
However, iHi for a (2£)-level system with only nearest neighbour interactions can have 
at most 2(2£ — 1) = 4£ — 2 non-zero entries on the first super- and sub-diagonal. Hence, 
a (2£) -level system with dynamical Lie algebra so(2£) must have interactions between 
no n- adjacent energy levels. 
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